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Long period Rayleigh wave and Love wave dispersion data, particularly for 
oceanic are^s, have not been simultaneously satisfied by an isotropic structure. 
In this paper available phase and group velocity data are inverted by a pro- 
cedure which includes the effects of traxisverse anisotropy, anelastic dispersion, 
sphericity, and gravity. The resulting models, for average Earth, average ocean, 
and oceanic regions divided according to the age of the ocean floor, are quite 
different from previous results which ignore the above effects. The models show 
a low-velocity zone with age dependent anisotropy and velocities higher than 
derived in previous surface wave studies. The correspondence between the 
anisotropy variation with age and a physical model based on flow aligned olivine 
is suggestive. For most of the Earth SH>SV In the vicinity of the low-velocity 
zone. Near the East Pacific Rise, however, SV>SH at depth, consistent witl« 
ascending flow. Anisotropy is as Important as temperature in causing radial andi 
lateral variations in velocity. The models have a high velocity nearly Isotropic 
layer at the top of the mantle that thickens with age. This layer defines the LID, 
or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum 
thickness of about 50 km. This thickness is comparable to the thickness of the 
elastic lithosphere, The LID thickness is thinner than derived using Isotropic or 
pseudo-isotropic procedures A new model for Average Earth is obtained which 
includes a thin LID, This model extends the fit of a P.R.E.M. type model to 
shorter period surface waves. 


INTRODUCTION 


Rayleig^ v^ave and Love wave phase and group velocity data have been col- 
lected and Inverted for upper-mantle structure by many authors. Some of the 
more recent studies are Schi'Ue and KnopofI [1977; 1978], Mitchell and Yu [1980], 
Silver and Jordan [1981], Forsyth [1975a; 1975b], Dzlewonskl and Anderson 
[1981] , Mills and Hales [1977], Wielandt and Knopoff [1982], Nakanishi and 
Arderson [1982], Anderson and Hart [1976], Nakanishi [1981], Anderson [1982; 
1983], and Montagner and Jobert [1983], Many studies have derived regional 
phase and/or group velocities for several regions believed to be physically 
different. For example, the earth can be divided into oceanic, tectonic, and con- 
tinental regions [Toksoz and Anderson, 1966; Nakanishi and Anderson, 1983; Mills 
and Hales, 1978b; Dziewonski and Steim, 1982; Silver and Jordan, 1981] or the 
ocean can be divided on the basis of the age of the seafloor [Forsyth, 1975a; 
1975b; Mitchell and Yu, 1980; Montagner and Jobert, 1983; Schlue and Knopoff, 
1977; 1978], However, the data derived in these studies have not been inverted 
using a procedure that accoiints completely and directly for the effects of aniso- 
tropy, as well as sphericity, gravity and anelasticity of the Earth. In this study 
an average Earth data-set, an average ocean data-set, and a data-set for a 
number of regionalized ocean age provinces will be modeled using a procedure 
that includes anisotropy, sphericity, gravity, and anelasticity. The data sets 
were collected from the literature. 

The regionalized oceanic data have been used to derive m.odels of the varia- 
tion with depth of SH velocity, using Love wave velocity data, and the variation 
with depth of SV velocity, from Rayleigh wave data. These models have been 
used to investigate how the velocity and the thickness of the oceanic lithosphere 
and the velocity in the low velocity zone (LVZ) vary with age. Forsyth [1975a] 
concluded that Rayleigh waves travel fastest in the direction of spreading and 
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that Love wave and Rayleigh wave data were inconsistent unless SH> SV in at 
least the upper 125 km. He also believed that anisotropy might be present 
below 250 km. Schlue and Knopofl! [1977] found that the LVZ was anisotropic but 
the crust and lid could be modeled as isotropic. Their LVZ e.xtends from 180 km 
to the bottom of the lid (15 km - 115 km depending on age). They comment that 
the observed velocity can be Cisplained by a thin sub-moho layer that would 
not be resolved by their data. Schlue and Knopofl [1978] Included anisotropy in 
the LVZ for a suite of calculated models. Yu and Mitchell [1979] and Mitchell and 
Yu [i960] find anisotropy predominantly in the lithosphere and possibly in the 
LVZ. These studies have found that lithosphere thickness and lithosphere velo- 
city generally increase with age, except Schlue and Knopofl! [1977] who con- 
strained their velocities to remain constant with age. The difierences in the 
depth where anisotropy is located in these studies can be attributed to 
differences in the assumptions, the constraints, the inversion methods, or to 
some systematic. 

The most serious source of systematic error is the separate isotropic inver- 
sion of Love and Rayleigh waves to give an anisotropic structure. The studies 
discussed above are almost all based on pseudo-isotropic inversions that deter- 
mine SV velocity from Isotropic Rayleigh wave inversion and SH velocity from 
Isotropic Love wave inversion. The differences between the two models are then 
used as a measure of the anisotropy. Separate isotropic inversions make no 
allowance for P-wave anisotropy, and include neither the effects of SV velocity on 
Love wave velocities nor the effect of changes in PV, PH, and SH velocities on 
Rayleigh waves. Thus, the procedure of using separate isotropic inversions is 
useful only to indicate the probable presence of anisotropy, or to calculate 
responses for propagation in planes perpendicular to the symmetry axis in a 
transversely isotropic medium [Crampin, 1976]. Some studies [Yu and Mitchell, 


1979; Moiitagner and Jobert, 1983j cite a statement made by Crampin and King 
[1977] that for low resolution data of an anisotropic earth an isotropic inversion 
is valid. Of course, the resulting model is necessarily of equally low resolution, 
and only the gross properties of the model may be believed. Thus the pseudo- 
isotropic inversion procedure is self limiting. No improvement in the amount or 
quality of the data will improve the quality of the model since the procedure 
itself becomes invalid when the resolution improves, Anderson [1966], Dziewon- 
ski and Anderson [1901], and Anderson and Dziewonski [1902] have shown that in 
the presence of anisotropy important errors are introduced by using a pseudo- 
isotropic procedure. Kirkwood [1970] also shows the necessity for the inclusion 
of anisotropy in the inversion. Thus, a complete anisotropic inversion would be 
useful and could produce different results than a pseudo-isotropic approach. 
An anisotropic inversion would have the advantagis that no error would be intro- 
duced by the incorrect pseudo-isotropic procedure, allowing increased resolu- 
tion in the presence of enough data. 

The inclusion of anisotropy introduces new difficulties to the problem. For 
the case of transverse anisotropy considered here the number of free parame- 
ters in each Layer increases from three (density, P velocity, S velocity) to six 
(density, horizontal and vertical P velocities, horizontal and vertical S velocities, 
and a velocity in some intermediate direction). The resolving power of the fun- 
damental mode surface wave data used here are insufficient to determine all of 
the parameters in models with even the simplest forms of pararrsterization. To 
overcome this problem additional constraints must be applied. The particular 
constraints used in this paper will be discussed later. 

As illustrated by the differences between the results of previous studies the 
parameterization, the class of models, and the constraints applied can have a 
significemt effect on the final models even for an isotropic inversion. While each 


of the studies was self-consistent and according to the appropriate set of resolv- 
ing kernels gave a good indication of the type and location oi anisotropy, the 
depth and extent of anisotropy was different for each study. While the 
difference In anisotropy may be partially caused by the pseudo-lsotroplc pro- 
cedure the differences between individual isotropic SV or SH models is also 
large. If one attributes those deferences to changes in model type, parameteri- 
zation, and constraints between studies, it is reasonable to expect that a full 
anisotropic parameterization, with its increased number of parameters, might 
yield different results. The increased number of parameters available vastly 
increases the number of possible models that will fit the data within a given 
range of uncertainty. The problem of non-uniqueness becomes even more pro- 
nounced, The final models are dependent on the parameterization and the 
assumptions. The models derived in this study are a possible set of anisotropic 
models that agree with the fundamental mode Rayleigh wave and Love wave data 
for periods less than 300 seconds. They can be used to illustrate how large a 
difference the introduction of anisotropy can make to models of the velocity 
structure of the upper mantle, 

DATA 

The data used in this study were compiled from many sources discussed in 
detail below. Data sets were collected for the average Earth, the average ocean, 
and for several regions of the Pacific divided according to the age of the ocean 
floor. Each data set consisted of phase and/or group velocities of Rayleigh 
and/or Love waves. For the average Earth, the average ocean, ana the regional- 
ized provinces 0-10, 0-20, 20-50, 50-100, and >100 M.Y. in age all four types of 
data were available. >’or the regionalized provinces 0-5, 5-10, and 10-20 M.Y, in 
age only Rayleigh wave data were available. 


The data for the average Earth are given in tables 1 and 2 and shown, along 
with the corresponding models, in figures I and 2, A large set of fundamental 
mode data was compiled from the sources discussed below. Nakanishi and 
Anderson [1983a; 1983b] give group velocities for Love and Rayleigh waves 
between 100 second.s and 400 seconds period and phase velocities between 100 
and 300 seconds period. Phase and group velocity values were determined from 
about 200 paths using the single station method or gf^eat -circle measurements. 
Spherically averaged velocities given here are part of the results of an inversion 
to give a spherical harmonic representation of the variations in great-circle 
velocity. Fukao and Kobayashl [1903] also give data of all four types for periods 
between 100 and 400 seconds. These data are the averages, and the quoted 
uncertainties are the standard deviy.tions of the velocities resulting from using a 
time variable filtering technique on data from 37 great-circle paths from the 
1963 Kurile islands earthquake. All fundamental mode data from these sources 
are included in the data set, Dziewonski and Steim [1982] give spheroidal mode 
periods and Rayleigh group velocities for qS’q to oS'gg that correspond to periods 
of 165 to 640 seconds. Selected modes for the period range 160 to 400 seconus 
are given in table 1. The modes chosen are q5i 7 to oS'ag, every second mode, and 
o^'ag to o'S'ssi every third mode. The intermediate modes are, for the most part, 
consistent with the other data and form a reasonably smooth curve with the 
modes used. Therefore, the included modes are a good indication of the quality 
of the fit of the model to the part of this data set with periods less than 400s. 
This data set was derived from 37 seismograms from several sources by using a 
waveform inversion technique to determine dispersion curves for each source- 
receiver pair. Then a pure path analysis was conducted to give group velocities 
and normal mode periods. Mills and Hales [1978a] give Rayleigh group velocities 
for periods of 50 to 300s, and normal mode periods for o'S'ss to o^’igg. Rayleigh 


group velocity values for periods between 50 and 300 seconds and phase veloci- 
ties derived from the normal mode periods of selected modes are given in table 
2. The uncertainties in the group velocity portion of their data set varied by an 
order of magnitude between two tables containing some of the same data. The 
uncertainties in Mills and Hales [1977] were taken as correct. The uncertainties 
quoted in Mills and Hales [1978a] were corrected to conform with Mills and Hales 
[1977]. These group velocities were derived ',-y averaging a set of path averaged 
great-circle group velocity estimates based on seismograms from several large 
events in the Kurile Islands. The phase velocities were determined by integrat- 
ing the group velocities. The data discussed above, although derived by different 
methods, are, in almost all, cases consistent between data sets to within quoted 
uncertainties. 

The data for the average ocean are given in tables 3 to 6, and shown in 
figures 3 and 4. Mitchell and Yu [1980] give data for four regionalized oceanic 
provinces. The data are derived from a pure path analysis of single station 
measurements for 33 Rayleigh and 30 Love wave paths in the Pacilic that are 
predominantly oceanic. To approximate average ocean the Rayleigh and Love 
group and phase velocity data for the two intermediate age provinces, 20-50 M.Y. 
and 50-100 M.Y., were averaged. This provides oceanic data for periods between 
20 and 110 seconds. Kanamori [1970] gives average oceanic Love and Rayleigh 
wave phase velocities between 125 and 300 seconds resulting from a pure path 
analysis of 25 great-circle phase velocity measurements. Mills .:nd Hales [1978b] 
give oceanic Rayleigh wave group velocities between 50 and 540 seconds. All 
average ocean data from these sources for periods less than 300 seconds are 
given in tables 3-6, and shown in figure 3 and 4. Dziewonski and Steim [1982] 
give oceanic Rayleigh wave group velocities and normal mode periods for aver- 
age ocean for periods between 165 and 635 seconds. These data were derived 
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using the analysis described in the average Earth section. The same modes used 
for the average earth data set are used for this data set. As for the average 
Earth data set, the remaining modes form a smooth curve with the modes used. 

Therefore, the modes used are a good indication of the quality of the fit to the 
model of the data set for periods less than 300 seconds. The collected data from 
the above sources are generally consistent with each other withi.i the quoted 
uncertainties. 

The data for regionalized oceanic provinces, predominantly from Mitchell 
and Yu [1980] and Forsyth [1975a; 1975b], are shown along with best fitting 
models in figures S to 8 and listed in tables 7 to 15, The regionalization of 
Mitchell and Yu [1900] consists of four provinces with age ranges 0-20 M.Y., 20-50 
M.Y., 50-100 M.Y., and >100 M.Y, For each of these regions short period data of 
all four types are given. Some of the regions of Forsyth [1975a, 1975b] were also 
used. These four regions are 0-5 M.Y,, 5-10 M.Y., 10-20 M.Y., and O-IO M.Y, The 
former data are derived from a pure path analysis of measurements in the 
entire Pacific basin while the latter data are derived from a similar study using 
measurements in the Nazca plate alone. For the 0-20 M.Y old region the East I 

Pacific Rise Rayleigh phase velocity data of Wielandt and Knopoff [1982] for 
periods of 40-300 seconds were added to the corresponding Mitchell and Yu 
[1980] data set that covers a period range of 20 to 102 seconds. The latter data 
.set is consistently slower by 0.005 to 0.02 km/s. The difference increases as the 
period becomes shorter. The differences are within quoted uncertainties of the 
data. Forsyth’s >20 M.Y. old province is equivalent in age to the 20-50 M.Y. old 
province of Mitchell and Yu [1980]. Combining the two sets of Rayleigh wave 
phase and group velocity data shows that in general Forsyth's data are about 
0.005-0.01 km/s faster. Again the differences are well within quoted uncertain- 
ties. The data sets for 50-100 M.Y. old ocean are entirely from Mitchell and Yu 
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[1900] and contain all four types of data. The data sets for >100 M.Y, oid ocean 
contained only Mitchell and Yu's data when the models were determined. How- 
ever, the Rayleigh phase velocity data of Souriau and Souriau [1983] were added 
after the model was derived. These data, which fit the model already caL'ulated 
to within two standard deviations, are based on a pure path analysis using the 
combined velocity data of Xanamori [1970], Dziewonski [1970], Wu [1972], 
Dziewonski and Gilbert [1972], Okal [1977], Nakanishi [1979], and Leveque 
[1980], Rayleigh wave data sets for the remaining provinces, 0-5 M.Y., 5-10 M.Y., 
10-20 M.Y., and 0-10 M.Y., are from Forsyth (1975a), Love wave data from For- 
syth (1975b) are given for the O-IO M.Y. old province. New Rayleigh wave group 
velocity data from Montagner and Jobert [1983] were compared with the models, 
and generally agree with these except in the 50-100 M.Y. old province. However 
.these data cover an increased period range of 40-300 seconds and thus would 
expand the data coverage, 

STARTING MODELS 

The introduction of anisotropy increases the number of free parameters in 
any inversion of the data and thus makes the choice of starting model and/or 
reference model more critical to the final model. The model parameterization 
and constraints need to be chosen to reduce the number of parameters the data 
will be required to resolve to a reasonable level. In each layer of a transversely 
isotropic model five elastic parameters, density, and a Q value, are needed to 
calculate a response. The elastic parameters can be in terms of the elastic con- 
stants, A, C, N, L, and F [Love, 1927] or in terms of wave velocities. The latter 
parameters were used for P.R.E.M. [Dziewonski and Anderson, 1901] and will be 
used here. Thus, the elastic parameters are VSV, VPV (vertical S-wave and P- 
wave velocity respectively), VSH, VPH (horizontal SH and P velocity respec- 
tively), and ETA. ETA determines how the P velocity and S velocity vary at 


intermediate angles of inuldence and can be expressed as follows: 
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To make comparison of the models to P.R.E,M. easier the structure of the 
P.R.E.M. model was used for this study, The models, therefore, consist of an 
ocean laver, an upper and lower crust, a region above the LVZ (the LID), a LVZ, 
and a region between the bottom of the LVZ (220 km) and 400 km. Below 400 km 
all models are Identical to P.R.E.M. As in the P.R.E.M. model anisotropy was 
introduced only between the base of the crust and the bottom of the LVZ; all 
other regions in the model are isotropic. The Q model from P.R.E.M. was used In 
all cases. The number of layers for which parameters must be determined was 
minimized by using a P.R.E.M. like polynomial representation. In the upper 400 
km, where parameters were allowed to vary in this study, each region (LID, LVZ, 
etc) can be treated as a single unit with a mean value and a linear gradient for 
each of the velocities, density, and ETA. Below 400 km P.R.E.M. has several 
regions described by linear, quadratic and cubic polynomials. These layers are 
not changed in any model discussed in this paper. Thus, the form of all models 
in this study can be described as follows: 

(1) :A water layer, a sediment layer, an upper crustal layer, and a lower crustal 

layer, all isotropic with constant velocities 

(2) : A LID whose thickness is allowed to vary at the expense of the LVZ, The LID 

has a thin (~6 km) upper layer with constant anisotropic velocities con- 
sistent with P-n observations. The remainder of the thickness of the LID is a 
layer of constant isotropic or anisotropic velocities and ETA, 
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(3) : A LVZ whose thickness decreases as the UD thickness increases. The LVZ Is a 

layer where the means and linear gradients of ETA and the four velocity 
parameters define the structure. 

(4) :A region between 220 and 400 kilometers where the linear gradients and 

means of ETA and the four velocity parameters define the structure. 

(5) : The density and Q in each region are the .same as In P.R.E.M. The means of 

the four velocities and ETA in all regions and gradients of velocities and ETA 
in regions 3 and 4 vary betv/een models. 

(6) : Below 400 kilometers the models are all identical to P.R.E.M. 

The starting models 'ised were all similar to P.R.E.M. but were different for 
each data set. Some features were common to many data sets. For example, 
the water and sediment depths for all the regional data sets were derived from 
those of Leeds, Knopoff, and Kausel [1974], A six kilometer thick upper Ud layer 
with velocities consistent with and Sn data was inserted below the crust for 
each of these regions, for the average Earth, and for the average ocean. Below 
400 km all models were identical to P.R.E.M. Above 400 km, density and Q in 
each depth interval were the same as in P.R.E.M. Differences in the starting 
models are discussed in the remainder of this section. 

For the average Earth data set P.R.E.M. itself was the starting model. For 
the average ocean the starting model consisted of P.R.E.M. below 40.8 km, a 30 
km LID with VPV=VPH=8.214 km/s and VSV=VSH=4.600 km/s and ETA=1.0, and 
the starting model of Yu and Mitchell [1979] for the crust sediment and water 
layers. 

For the regionali 2 .ation of Mitchell and Yu [1980] starting models were also 
similar to P.R.E.M. For 0-20 M.Y. old ocean a 3.45 i. n water layer and a 0.05 km 
sediment layer were placed on the same crustal starting model as the average 
ocean data. A 14 km thick lower lid layer with VSV=ViTN."4.6 km/s and 


VPV=VPH=8.2 km/s was inserted and the LVZ velocities were reduced by 3.5% 
with respect to P.R.E.M. For the 20*50 M.Y, old ocean a 4.4 km water layer and a 

O. 1 km sediment layer were placed over the same crust, a 19 km lower lid layer 
with ySV=VSH=4.68 km/s and VPV=VPH=8.35 km/s, and a LVZ with velocities 
reduced 2.7% with respect to P.R.E.M. For the 50-100 M.Y. old ocean a water 
layer of 4,7 km and a sedimevH layer of 0.3 km as given by Yu and Mitchell 
[1979] were used. Latoi--, the water layer thickness was adjusted to 5,4 km and 
the sediment to 0.23 km in a modifled starting model consistent with Leeds, 
Knopoff and Kausel [1974]. A 34 km thick lower lid with VSV=VSH=4.B km/s and 
VPV=VPH=8,57 km/s, a LVZ with velocities decreased by 2.25% with respect to 

P. R.E.M., and a region between 220-400 km were velocities increased by 1% were 

used. For ocean >100 M.Y, in age a water layer 5.75 km thick, a sediment layer 
0.3 km thick, a lower lid 39 km thick with VSH=VSV=4.75 km/s .-ijnd a LVZ with 
velocities increased by 1% and ETA increased by 1.2% wiib to P.R.E.M. 

For Forsyth's regions the final modc4 for the 0-20 M.Y, old ocean (table 16) 
with minor change was used as a starting model, The lid thickness was set to 9 
km for 0-5 M.Y., 17 km for 5-10 M.Y., 14 km for 10-20 M.Y. and 0-10 M.Y. The 
water depth was set to 3.3 km for 0-5 M.Y., 3.5 km. for 5-10 M.Y., 3.8 km for 10-20 
M.Y,, and 4 km for 0-10 M.Y.. No sediment layer was included except for 10-20 
.M.Y. where 0,02 km of sediment was added. 

PROCEDURE 

The velocity depth models discussed in later sectionfs were determined by 
iterative forward modeling. The procedure used to calculate the dispersion 
curves treats transverse anisotropy in a spherical, anelastic, self gravitating 
Earth [Takeuchi and Saito, 1972; Dziewonski and Anderson, 1981]. Periods of 
spheroidal and torsional modes of free oscillation are determined for the given 
structure. The phase velocities for Rayleigh and Love waves are then calculated 
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from. 


r X (71 + 

where T is the period, n is the mode number, oS'n or oTr. and is the radius of 
the earth. For all calculations we use = 6371 km, 

Modes were chosen to minimize computation time while supplying points in 
the velocity-period plane separated by not more than twenty .seconds in period. 

0 

The modes calculated for each iteration were; o-^as to o‘S' 40 i every third mode, 
to o5g8. every tenth mode, o54e to qS^zq. every thirtieth mode, qTzq to o7’4e< 
every third mode. o7’4g to oFios. every tenth mode, and o^ios to oT 4 gs, every thir- 
tieth mode. Linear interpolation between, these points gives the velocity values 
at the periods corresponding to data points, The final models were calculated 
for more densely packed modes so Interpolation to the data would be more 
accurate. 

Changes in the models are chosen to reduce the deviation between data and 
model. At each iteration at least one and usually not more that three p,arame- 
ters were changed. Occasionally several parameters were changed simultane- 
ously by the same relative amount. This has the effect of varying one parameter 
and requiring that its ratio with each of the other changed paranaeters remain 
constant. 

Elach of the data sets was fitted by making minimal changes from the 
appropriate P.R.E.M. type starting model. An attempt was made to keep the 
variations of parameters wi.th age smooth but as an iterative forward modeling 
technique was used no formal smoothness criteria was applied. To keep the 
number of parameters small, an experiment where parameters were succes- 
sively freed was conducted and is described in the following paragraphs. The 
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efl’ect of each of the possible parameters on the dispersion curves was calcu- 
lated and the tradeoffs between parameters noted. This will be discussed later, 
Considering the number of parameters, the correlations between them, and the 
limited resolution of the data, the non-uniqueness resulting from the under- 
determined nature of the problem will make physical Interpretations of the 
resulting models somewhat speculative. However, this type of procedure can be 
useful in testing the validity of any theoretical upper mantle model that can be 
suitably formulated. 

To determine a proper parameterization a series of experiments with 
different parameters was performed. The average ocean data of Mitchell and Yu 
[1980] were modeled by successively freeing each parameter. LID thickness and 
isotropic LID velocity were varied to produce a set of parametric curves in the 
two parameters. A six kilometer anisotropic layer was placed at the top of the 
lid to make the model consistent with observed velocities and to improve the 
fit at short periods. To reduce model velocities to values near the data the mean 
SH velocity in the LVZ was decreased while VSH/VSV, VSH/VPH, and VPH/VPV 
were held constant, This magnified the changes due to lid thickness and velo- 
city. The magnification is caused by the increased contrast in velocity across the 
LID LVZ discontinuity. To resolve incoL.nstencies which prevented simultaneous 
fitting of Love and Rayleigh wave datei, particularly at periods below 100 seconds, 
the velocities in the LID and the LVZ were decoupled such that only VSV/VPV and 
VSH/VPH ratios remained constant, that is VSV/VSH was allowed to vary. Then 
the gradients of velocity and ETA were varied in the LVZ and the means and gra- 
dients of velocity were completely decoupled. To improve the fits at shorter 
period.s the crustal velocity, crustal thickness, water thic.kness, and sediment 
thickness were varied. These variations were small ,a couple of percent max- 
imum, but they dramatically improved the fit for periods less than about 80 
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seconds. Thus the parameters investigated are these: 

(t); water depth 

(2) : sediment thickness 

(3) : crustal velocity: the crust is assumed to be isotropic and the velocities 

in the two layers are varied by the same relative amount. 

(4-5); crustai thickness; the two crustal layers are separately varied 
(6): lid thickness; 6 km upper layer , thickness of remaining lid to be deter- 

mined 

(7-8): lid velocity: VSV/VPV and VSH/VPH are held constant, VSH and VSV 

are allowed to vary in the lower lid layer, velocities in the upper lid 

layer are fixed to be consistent with data. 

(9-13); mean VSH, VSV. VPH, VPV, and ETA in the LVZ. ' 

(14-18): gradient of VSH, VSV, VPH, VPV in the LVZ 

(19): Mean VSH between 220 and 400 km. VSH APH, VSV/VPV, and VSH/VPV 

are held constant. 

Each of the nineteen parameters discussed above produce specific changes 
to the phase and group velocity dispersion curves. There exist several tradeoffs 
between these parameters that make the models non-unLque within the error of 
the data. However, some tradeoffs can be limited by the application of con- 
straints. The partial derivatives of the dispersion curves with respect to each of 
the above parameters were determined by differencing the Love and Rayleigh 
phase and group velocities calculated in two almost identical models. The two 
models differed in only one parameter. The resulting estimates of the partial 
derivatives are shown in figures 9 to 12. The correlations between them and 
their behavior are discussed in detail below. 

Variations in parameters associated with the water layer, the sediment 
layer, and the crustal layers arc interrelated. Varying the depth of the water 




and S(vdiment layers causes the curvatiUw* 'if tht dispersion curves to change at 
periods less then 45 seconds. For an increase in the water or sediment thick- 
ness the reduction in phase velocity increases as period decreases since short 
periods are more sensitive to the near surface velocity perturbations. Crustal 
thickness was increased by about 6.5%, and crustal velocity was increased by 5%, 
Increasing the crustal thickness at the expense of the faster LVZ means a net 
decrease in velocity in the depth range between the bottom of the old crust and 
the top of the new LVZ. Thus, the sharp decrease in group and phase velocities 
at periods less than 45 seconds is reasonable as is the sharp increase in the 
same period range when the velocity of the crust is increased. The partial 
derivatives in figures 9 to 12 show that comparable effects at short period result 
when crustal velocity increases and when crustal thickness decreases, although, 
the former has a larger effect at periods of 50-100 seconds. Also a similar effect 
occurs at short period(<45 seconds) when water depth is decreased. Observed 
values of water depth and crustal thickness and velocity are used to bound the 
models derived in this study. 

The effects of variations in the parameters associated with the lid are 
changes in the dispersion curves that are most obvious at periods less than 100 
seconds but are easily measurable at all periods considered. Increasing VSV and 
VPV in the LID causes a small but sharp increase in Love wave phase and group 
velocity at periods <60 seconds and a steep increase in Rayleigh velocities for 
periods <100 seconds. Decreasing VSH and VPH in the LID causes a small 
decrease in Rayleigh velocities and a larger decrease in Love velocities. It is 
important to note that changes in LID velocity must be weighted by the LID 
thickness to determine the correct perturbations to improve the fit of the 
model. The partial derivatives shown are for a increase in the mean velocities of 
3% for SV and PV and for a decrease in the mean velocities of 5% for SH and PH. 
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The effects of variation of the five parameters associated with the mean 
velocities and mean ETA in the LVZ are strongly interrelated. The mean of each 
velocity in the LVZ was increased about 5% and the mean of ETA was increased 
about 2%. This translates to an increase of 0.4 km/s in P velocity, 0.2 km/s in S 
velocity, and 0.02 in ETA. The resulting partial derivatives show that Rayleigh 
waves are most sensitive to VSV, as expected, but they also show that Rayleigh 
waves are sensitive to VPV and VPH and ETA at a significant level, while the effect 
of VSH is an order of magnitude smaller. In contrast the Love waves are most 
sensitive to VSH, as expected, but are also sensitive to VSV. This illustrates one 
of the problems of the pseudo- isotropic procedure. In an anisotropic earth Ray- 
leigh waves have a significant dependence on P velocities and on VSH velocity 
and Love waves have a significant dependence on SV velocity. The problem is 
coupled rather than separable as required by the pseudo-isotropic procedure. 

The effects of variations in the five gradient parameters associated with the 
LVZ are strongly interrelated and are intimately related to the effects for 
changes in the mean parameters. The effects of changing the gradients in the 
LVZ while keeping the means constant was to produce partial derivatives of the 
same sign as changing the corresponding mean for periods greater than about 
100 seconds, and to produce partial derivatives of the opposite sign for periods 
less than 100 seconds. The relative magnitudes of the contributions of each of 
the velocity gradients was in proportion to the contributions from the 
corresponding means. Tradeoffs exist between groups of velocities and velocity 
gradients in the LVZ but they involve complex combinations of several parame- 
ters. To illustrate the possible tradeoffs a discussion of two equally well fitting 
models to the average oceanic data is presented below. Decreasing the thick- 
ness of the LID will increase the thickness of the LVZ and will proportionately 
increase the effects of the parameters in the LVZ. This means that the 
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thickness of the LID has more effect than the direct partial derivative in thick- 
ness indicates., 

The procedure used to forward model each data set considered magnitudes, 
shapes, and interrelations of the partial derivatives of the parameters. In order 
to find a model to fit each data set within a specified uncertainty while adjusting 
the smallest number of parameters the following procedure was used: 

(1) : Assign values for water, sediment and crustal parameters appropriate for a 

given data set (eg. water and sediment thicknesses from Leeds KnopofI and 
Kausel (1974), or crustal structures from Yu and Mitchell (1979), and For- 
syth (1975a)) 

(2) : If necessary adjust mean SH velocity in the LVZ or between 220 and 400 km, 

keeping VSH/VPH, VSH/VSV, and VSV/VPV ratios constant, to fit the longer 
periods. To keep the variation of velocities with age as smooth as possible 
the magnitudes of changes in the two regions can be traded off, 

(3) : Adjust the thickness, the velocities, and ETA of the lid to fit the data as well 

■is possible. 

(4) : If necessary adjust the gradients of velocity and ETA in the LVZ keeping all 

ratios constant. 

(5) ; If necessary decouple the means and gradients of the velocities and ETA in 

the LVZ such that VSH/^PH and VSV/d/PV remain constant and VSH/VSV is 
allowed to vary. 

(6) : If necessary decouple VSH/VPH and VSV/VPV 

(7) : If necessary fine tune by returning to 3 and repeating the process until the 

fit is at the desired level of accuracy. 



(8): Fine tune the fit at short periods by varying water, sediment, and crustal 
velocity. The changes should be small if any are needed. 

The problem of non-uniqueness is particularly obvious v^hen using this type 
of method. Very different models can fit the data to within comparable standard 
deviations. As an example of this, two models that fit the average ocean data 
with comparable standard deviations are discussed below. The two velocity 
depth profiles are shown In figure 13. Except for the SH velocities the models 
are quite different. There are substantial differences in the means of the other 
three velocities in the LVZ and the corresponding three velocity gradients differ 
in sign. The gradients of ETA are also different. While the fits of these two 
models to the data, shown in figures 3 and 4, are not Identical the standard devi- 
ations of the two sets of residuals are comparable. The Rayleigh residuals are 
given in table 3 and 4- and shown if figure 14. The hove residuals are given in 
tables 3 and 6 and shown in figure 15. 

The standard deviation of the residuals are given in table 18; S is the stan- 
dard deviation giving each point equal weight, SW is a weighted standard devia- 
tion where the weights are proportional to the distances in period between adja- 
cent points. The latter statistic should remove any bias caused by unequal dis- 
tribution of data points in period, and is determined from the following rela- 
tions: 


Sfy^ = 


tnn - 

f \ 

n 

2 / 1 


i=\ 

[i=l 



S/i -1 

i = \ 



where 


- 20 - 


OF 0'^ v.Ii'V 


/i = 


ra-r, 


/n ” T'n“"‘'^n - I 


f l+i = 


7’i.i-T’i-i 

^^2 


for -i = l, . . . ,n-l 


For Rayleigh phase velocities model B fits better at shorter periods and 
model A fits better at longer periods. Overall, model B is a slightly better fit but 
the difference of about 1.4% should not be significant. For Rayleigh wave group 
velocities model A gives a significantly better iU. For Love wave phase velocities 
model B gives a marginally better fit and for Love wave group velocities the qual- 
ity of fit is essentially identical. A small change in the SH velocities of model A 
would increase the Love residuals so that the two model would give equal stan- 
dard deviations, and such a change increases the difference between the two 
models. Thus, quite different models can fit the data to the same accuracy, 

DISCUSSION 

The models derived using the method described above are shown in figures 
16 through 20. Calculated Rayleigh and Love wave phase and group velocities 
are listed in tables 1 through 14. In this section the models resulting from each 
of the data sets and some possible implications arising from them will be dis- 
cussed, The most striking difference from previous studies, common to all 
models, is the thickness of the lithosphere. The lithosphere , defined here as the 
LID plus the crust, is much thinner than the lithosphere derived in. isotropic 
seismic studies and is comparable in thickness to the flexural lithosphere 
defined by Watts et. al. [1900], Also, the variation of anisotropy with age in the 
regionalized oceanic data sets indicate that this method can be used to test 
models of mantle flow based on the hypothesis of aligned olivine crystals. The 
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regionalized oceanic data sets of \fitchell and Vu [1980] and Forsyth [1975a] 
both show an increase in lithosphere thickness with age, and variations in litho- 
sphere and LVZ velocities and ETA values with age. 

The average Earth velocity model is presented in figure 18 and in table 17, 
the phase and group velocities at the periods of the data are given In tables 1 
and 2. The average Earth model is similar to P.R.E.M. Above 18.4 km and below 
220 km no changes have been made. The major difference is the introduction of 
a constant velocity LID region between 18.4 and 46.8 km depth. This causes a 
large increase in velocity between the bottom of the P.R.E.M crust at 24 km and 
the bottom of the average Earth model crust at 18.4 km. Thus, the mean veloci- 
ties and ETA values over the depth range of the LID are Increased by about 2-4% 
for velocities and decreased about 2.7% for ETA. The introduction of the LID 
reduces the depth extent of the LVZ. Mean velocities in the LVZ increased by 
0. 1-1.5% for SH, PH, and PV and decreased by 0.25% for SV, The range of the 
velocities and ETA values in the LVZ Increase by about 12% for SH and PH to 
0.018 km/s and 0.015 km/s respectively and decreased considerably for SV and 
PV, to about 0,012 km/s. These changes provide a better fit than P.R.E.M. to 
the short period fundamental mode data. 

Figures 1 and 2 show that this model fits the fundamental mode data well. 
Overall the residuals have a mean of -0.0031 i 0.010'2 km/s. The residuals, given 
in tables 1 and 2, show that data from different sources have systematic 
differences. The model consistently lies below the data of Fukao and Kobayashi 
[1983] except for shorter period Love phase velocities. It is within the quoted 
uncertainty of the data 96% of the time. The data of Nakanishi and Anderson 
[1983] were added after the model was derived. They show residuals that are 
generally smaller than the other data sets. However, the extremely small stan- 
dard deviations of the data cause more than half the points to be separated 
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from the model by more than 2a, The model lies consistently above the Rayleigh 
phase velocity data of Mills and Hales [1977]. Their Rayleigh wave data set lies 
within their quoted uncertainties 33% of the time. 

The average ocean model is shown in figure 13 and given In table 15, The 
model shown in dotted lines, model B, is most similar to P.R.E.M, The dispersion 
curves resulting from this model are shown in figure 4, and the phase and group 
velocities at the periods of the data points are given in tables 3 to 6. There are 
several major differences between P.R.E.M. and the model for average ocean. 
First, the water layer is thicker since it is not being averaged over the entire 
surface of the Earth, Also, the crustal layers are thinner since the thicker con- 
tinental crust is not being Included in the average. Therefore, there is a net 
decrease in velocity of 13.6% for S and 9.1% for P caused by the additional thick- 
ness of the slow water layer and to a lesser extent to the decreased thickness of 
the faster crustal layers. Finally, a LID was introduced, between 12 km and 53 
km. Since, between 12 km and 24.4 km the average ocean LID corresponds to 
the slower crustal layers in P.R.E.M. the oceanic model's LID Is faster by 11,4% 
for PV, 10.1% for PH, 12.0% for SV, 9.3% for SH, and 5,0% for ETA when compared 
to the mean velocities over the same depth ranges in P.R.E.M. In the LVZ the 
means of velocity and ETA increase slightly. For PV, PH, and SH slopes became 
steeper, and for SV and ETA slopes were reduced. 

Figure 4 shows how this model fits the fun dt- mental mode data set. Again 
the different data sets show systematic differences. The data of Mitchell and Yu 
[i960] are well fit with the Love wave velocities and the Rayleigh wave phase 
velocities lying consistently above the data and the Rayleigh wave group veloci- 
ties lying consistently below the data. The data of Kanamori are well fit by the 
model. All points lie within 1.1 a. The Rayleigh wave group velocity data of 
Dziewonski and Steim [1983] all are fit by the model at a 2a level, and 55% lie 



within la. The model lies consistently above the Rayleigh wave phase velocity 
data of Mills and Hales [1970a] and fits at the 3a level. Overall, the residuals 
have a mean of 0.0024±0.0312 km/s, with 81% of points failing within la. Thus, 
this is a good fit to the data. 

The models for the regions used by Mitchell and Yu [1980] and Forsyth 
[1975a] show an Increase in the thickness of the lithosphere with age. Particu- 
larly for the younger provinces of Mitchell and Yu [1900] the models fit much 
better than those previously proposed. The models were derived by trying to 
make a smooth variation with age in as many parameters as possible, The aver- 
age ocean model included in this series of models Is model A which fits smoothly 
into the progression of parameters in the regional models. This average ocean 
model was derived using only the data of Mitchell and Yu [1900]. It is reasonable 
to assume that a set of models similar to the alternate average ocean model 
could be derived to fit as well as, or better than, the models presented here. 

The velocity depth models for the 0-20 M.Y. and 20-50 M.Y. old. oceanic 
regions are shown in figure 17 and given in table 15, the phase and group veloci- 
ties and the residuals are given in tables 7 and 8, and the corresponding disper- 
sion curves are shown in figure 5. The data of Mitchell and Yu [1980] are well fit. 
For the younger region Love group velocities tend to fall below the model at 
short periods and above it at longer (>70s) periods. The model lies consistently 
1-2 a below the Rayleigh wave phase velocity data of Wielandt and KnopofT for 
periods >180s. Below that period the data is fit at a la level and there is no sys- 
tematic bias. For the older region the model fits both data sets well. 

The velocity depth models for 50-100 M.Y, and >100 M.Y. old oceanic regions 
are given in figure 10 and table 15, corresponding velocities are given in tables 9 
and 10, and dispersion curves are shown in figure 6, For both data sets the 
phase velocities are very well fit and the model tends to lie above the data for 


longer period Iiove phase velocities and below it for shorter periods. Group velo- 
cities do not flt as v'ell but still agree well with the data at a 2a level, The data of 
Souriau and Souriau [1983] added after the model was derived, fit at a 2a level 
and add data coverage between 102 and 300 seconds for the >100 M.Y, old 
region. The uncertainties of group velocity points at a given period are smaller 
for the older regions. The poorer fits mentioned above are judged with respect 
to those uncertainties. Thus, although the absolute residuals are comparable 
the agreement between data points and the model appears to be poorer for 
older regions where uncertainties are smaller. 

The velocity depth models for 0-5 M.Y. and 5-10 M.Y, vld oceanic regions are 
given in figure 19 and table 18, the corresponding dispersion curves are shown in 
figure 7, and the velocities and residuals are given in tables 11 and 12. The velo- 
Q’ly depth models for 10-20 M.Y. and 0-10 M.Y old oceanic regions are given in 
figure 20 and table 16, the corresponding dispersion curves are given in figure 8, 
and the velocities and residuals are given in tables 13 and 14. These four models 
flt the data quite well. There appear to be no systematic trends in the relation- 
ship of data and models. 

The differences between the regionalized models can be interpreted as a 
progression of some physical properties with age. First, the thickness of the LID 
increases with age from 16 km for 0-5 M.Y. old oceanic regions of the Nazca plate 
to 26 km for the corresponding 10-20 M.Y. old regions, and from 20 km in the 0- 
20 M.Y. old oceanic regions of the Pacific to 50 km in the Pacific oceanic region 
>100 M.Y. in age. A somewhat thicker LID might be accommodated by adjusting 
the other parameters but inserting a LID of the thickness found in most other 
studies would be difficult, if possible at all, using this approach. The differences 
in lithospheric thickness between typical previously accepted models and the 
models presented here are small for young regions and increasingly serious as 


the ocean floor increases in age. The variation of LID thickness with age is shown 
in figure 21. This figure also shows the estimates by Watts et. al. [1980] of the 
thickness of the -elastic or flexural lithosphere. Next, the lid velocity increases 
with age for all but the two oldest provinces. In these two oldest regions the lid 
velocity becomes anisotropic with vertical velocities slower than horizontal velo* 
cities, Also, the water and sediment layers increase in depth with age, the cru- 
stal thickness is allowed to increase slightly after 20 M.Y., and the crustal veloci- 
ties decrease slightly after 50 M.Y Up to 20 M.Y. age the crustal thicknesses 
and velocities are unchanged. Thereafter, changes are small and result from 
fine tuning the fits at very short periods (<50s), The dispersion data, for old 
ocean, could alternately be fit by further increasing crustal and sediment thick- 
ness, or water depth. However, by using reasonable observational estimates as 
constraints on these thicknesses the velocity decrease in the crust seems neces- 
sary, The gradients of SH and PV velocity increase with age. .Many of the trends 
such as increasing lid velocity, increasing SH and PH velocities at the Lop of the 
lid, etc. are not continued into the >100 M.Y old province of Mitchell and Yu 
[1930]. A possible reason for this is the presence of many oceanic plateau areas 
in the old ocean region. Inclusion of such slow structures would reduce the 
average velocities for the region causing the progre-s-sion of velocities with age to 
be broken. The lateral heterogeneity introduced by plateau regions might also 
account for the need to introduce anisotropy in the lid of the homogeneous but 
anisotropic model. 

The thickness of the lithosphere is an important parameter in the theory of 
plate tectonics. The seismic lithosphere, sometimes called the LID, has previ- 
ously been calculated to be at least 100 km thick for oceanic regions >100 M.Y. 
in age. This value is commonly accepted and is often adopted as the thickness 
of the slab, or plate. The plate is considered to translate coherently in most 
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plate tectonic models, Thermal cooling models of the aging plate have also 

sometimes been based on this thickness. In turn these thermal models are used 

to determine thermal properties of the LID or, if thermal properties are 

assumed, to determine the thlekness of a thermal lithosphere. The thickness of 

the plate is also Important when r.aodellng convection and other physical 

processes in a subduction zone. Kowever, in subduction zones there are some 

independent evidences for thinner lithosphere, For example a thin LID was 

inferred from long range refraction measurements in the western Pacific 

[Nagumo ct. al,, 1981; see triangle in figure 21], It is also interesting to note 

that the separation of the earthquakes in the double Wadati-Benioff zone in 

Japan [Hasegawa et. al., 1978], if interpreted as the thickness of the subducting 

LID, is comparable to the LID tfl‘'ckness derived in this study. It has also been 

observed that the lithosphere in the region of some hot spots is thinner than 

predicted by the age of the lithosphere in question. This has been interpreted 

as a thermal resetting of thickness due to the hot spot [Detrick and Crough, 

1978], However, these measurements could also be used to support the 

hypothesis that the seismic lithosphere is thinner than previously believed, as 

derived by anisotropic modelling. If further studies of other types of data can 

verify the thicknesses found for the models presented here, and by Anderson 

and Regan [1983], then many accepted interpretations and ideas will need to be 

reconsidered. 

* 

The thickness of the elastic lithosphere has been determined by studying 
phenomena such as seamount loading [Watts et. al., 1980]. Generally it is 
believed that the elastic lithosphere rapidly thickens with age to a thickness of 
about 30 km for seafloor 50 M.Y. in age. Thereafter, the thickness seems to 
increase more slowly. Previous estimates of the thickness of the seismic litho- 
sphere, based on isotropic calculations, agreed that it was considerably thicker 
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than the elastic lithosphere, The difference was reconciled by discussing the 
different phenomena involved in determining the thicknesses. Each phenomena 
was discussed in terms of reiaxation times [Anderson and Minster, 1980] for 
relaxation by dislocation glide and climb. For this type of interpretation the 
thickness of the iithosphere is defined as the depth having a characteristic time 
equal to the duration of the load. Relaxation times generally decrease with 
depth as temperature increases. For high stress, long duration loads, .such as 
seamount loading and post glacial rebound, the thickness of the lithosphere is 
small and the relaxation time is long. For low str'iss short duration seismic 
waves relaxation times in and and for some distance below the thin flexural 
lithosphere are long compared to seismic periods. Thus, the seismic lithosphere 
is thicker than the elastic lithosphere in a homogeneous mantle. The relaxation 
time is strongly dependent on the temperature and stress. Thus, the thickness 
of the flexural or seismic litnosphere is highly dependent on the temperature 
and stress profiles. 

The interpretation discussed above is based on dislocation glide and climb 
which are both thermally activated processes. Thus, the relation between 
flexural and elastic lithosphere thicknesses is dependent on temperature. How- 
ever, the strength and other rheological properties of the lithosphere also 
depend on mineralogy, crystal structure and orientation, stress, partial melting, 
and duration of load. If temperature is not the dominant parameter, as could be 
the case if the upper mantle is chemically layered with the base of the LID 
corresponding to a change in crystal structure or mineralogy, the elastic and 
Seismic lithosLpheres might be of the same thickness. 

The variation of velocities and anisotropy with age suggests that interpret- 
ing the results in terms of stress or flow aligned olivine [Nicolas and Poirier, 
1976] is a promising approach [Regan and Anderson, 1981]. Using 'such an 


approach the velocities depend primarily on temperature, pressure and crystal 
orientation. The effects of each of these variables will be discussed below. The 
basis of an interpretation based on flow and the resulting velocity depth rela- 
tions are illustrated in flgure 22. The upper left diagram schematically illus- 
trates a convection cell with material rising at the midocean ridge (R) and 
flawing down at the trench (T). 

In the lower left of flgure 22 is the schematic temperature profile for such a 
cell. The seismic velocities decrease with temperature and increase with pres- 
sure as depth increases. Combining the effects of temperature and pressure 
one obtains a relation between velocity and depth. At the ridge the temperature 
increases very rapidly with depth near the surface. Thus, the effects of tempera- 
ture dominate over those of pressure and velocities decrease. Deeper levels 
under the ridge are almost l.sothermal. Thus, the effect of pressure dominates 
and the velocities increase. At the trench the temperature gradient is large 
near the base of the cell and nearly isothermal at shallower depths. Therefore, 
the velocity response is a mirror image of that at the ridge. Midway between the 
ridge and the trench (M) the temperature increases rapidly near the top and 
bottom of the cell. Thus, the velocities decrease rapidly in these regions. 

The crystal orientation, if alignment with flow is assumed, is with the shor- 
test and slowest axis (b-axis) perpendicular to the flow. Thus, at the ridge and 

I 

the trench where the flow is vertical the b-axis is horizontal, and midway 
between where the flow is horizontal the b-axis is vertical. The following discus- 
sion considers the effects of crystal orientation alone in a flowing olivine aggr'^- 
gate. VPH and VSH midway between the upward and downward flowing edges of 
the convection cell are controlled by the velocities along the a-axis and c-axis; 
VPV and VSV depend on the velocity along the b-axis . Thus, at the midpoint, 
and wherever flow is horizontal, SH>SV and PH>PV. At the ridge and at the 


trench the flow is vertical, rapidly changing to horizontal at the top and bottom 
of the cell. For vertical flow the horizontal velocity Is controlled by the b-axis 
and c-axis velocities so SH<SV and PH<PV. The values at the top and bottom of 
the cell rapidly change to the horizontal flow values, Between the midpoint and 
the trench or ridge the transition from horizontal to vertical flow velocities 
becomes sharper and the depth extent of constant vertical velocities increases. 
Combining the effects of temperature and pressure with the effect of orien- 
tation requires a calculation of the magnitude of each of the effects. The values 
of the elastic parameters and their derivatives with pressure and temperature 
have been measured for olivine and fosterite [Graham and Barsch, 1969; Kuma- 
zawa and Anderson, 1969], and for olivine rich rocks [Christensen and Crossen, 
1968; Christensen and Smewing, 1901], To determine the size of the velocity 
gradients caused by pressure and temperature effects these measured values 
were substituted into the following relation that defines the temperature deriva- 
tive of velocity, 

dV _ dV dP dV dT 

dT ~ dP dz dT dz 

9 7 ^ 

The temperature profile shown in fig 22 was used to determine ~ — : The meas- 

dz 

ured values of the elastic constants give the magnitude of the orientation effect. 
Thus, the two effects can be combined to produce velocity profiles like those 
shown in figure 22, 

Comparing this model to the derived velocity depth structures for the 
regionalized oceanic provinces shows that it could be a viable interpretation. 
For the average ocean model and for the upper 100 km of the mantle in the 
youngest regions PH>PV and SH>SV, This is consistent with horizontal flow. For 
the 100-220 km depth range for the youngest regions (0-5, 5-10) SH>SV The 
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vertical flow expected in the ridge crest environment would exhibit this 

dV 

behavior. The temperature gradients implied by the measured values of ^^and 
dV dV 

and the values of — - from the models are 5-8 degrees Centigrade per 
0 1 dT 

kilometer for older ocean and are consistent with reorientation of olivine along 
with a small temperature gradient for younger oceanic regions. Note that with 
these temperature and flow models the velocity of Love waves along ridges is 
extremely slow and tbe velocity of Rayleigh waves is high along subduction 
‘/lones, For midplate locations Love wave velocities are higher and Rayleigh wave 
velocities are lower than at plate boundaries. The results of this study Indicate 
that a model derived on the basis of flow aligned olivine could be fitted by exist- 
ing data with small changes to the present models. 
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200,00 

4.5880 

4.5870 

4.5735 

0.0130 

-0.0010 

-0.0145 

3 

190.34 

4.5228 

4.5268 

4.5127 

0.0041 

0.0040 

-0,0101 

3 

180.99 

4.4681 

4.4734 

4.4584 

0,0044 

0.0053 

-0.0097 

1 

175.00 

4.4420 

4.4391 

4.4236 

0.0110 

-0.0029 

-0.0184 

3 

172.54 

4.4194 

4.4251 

4.4093 

0,0048 

0.0057 

-0.0101 

3 

184.83 

4.3758 

4.3811 

4.3645 

0.0057 

0.0053 

-0.0113 

1 

150.00 

4.3250 

4.3045 

4.2871 

0.0320 

-0.0205 

-0.0379 

2 

102.40 

4.0820 

4.0994 

4.0815 

0.0180 

0.0174 

-0.0005 

2 

93.10 

4.0510 

4.0690 

4.0521 

0.0190 

0.0180 

0.0011 

2 

81.90 

4.0250 

4.0406 

4.0258 

0,0190 

0.0156 

0.0008 

2 

70,60 

4.0050 

4.0150 

4,0026 

0.0190 

0.0100 

-0.0024 

2 

60.00 

3.9940 

3.9996 

3.9900 

0.0190 

0.0056 

-0.0040 

2 

51.20 

3.9910 

3,9940 

3.9870 

0.0190 

0,0030 

-0.0040 

2 

40.00 

3.9960 

3.9967 

3.9927 

0.0200 

0.0007 

-0.0033 

2 

35.30 

4.0010 

4.0010 

3.9978 

0,0200 

0,0000 

-0,0032 

2 

60.00 

3,9940 

3.9996 

3.9900 

0.0190 

0.0056 

-0.0040 

2 

51.20 

3.9910 

3.9940 

3.9870 

0.0190 

0.0030 

-0.0040 

2 

40.00 

3.9960 

3.9967 

3.9927 

0,0200 

0,0007 

-0,0033 

2 

35.30 

4.0010 

4.0010 

3.9978 

0.0200 

0.0000 

-0.0032 

2 

30.00 

4.0060 

4.0073 

4.0044 

0.0200 

0.0013 

-0.0016 

2 

28,00 

4.0070 

4,0094 

4.0063 

0,0200 

0,0024 

-0.0007 

2 

26.00 

4.0060 

4.0106 

4.0073 

0.0200 

0.0046 

0.0013 

2 

24.00 

4,0040 

4.0099 

4.0063 

0,0200 

0.0059 

0.0023 

2 

22.00 

3.9970 

4.0051 

4.0014 

0,0200 

0.0081 

0.0044 


References: 

N=1 Kanamori (1970) 

N=2 Mitchell and Yu (1980) 
N=3 Mills and Hales (1978b) 
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! RAYLEIGH WAVE GROUP VELOCr 

PABLE 4 
PIES FOR ' 

PWO AVERAGE OCEAN MODELS 

T" 

T 

Co 

Cb 

Ca 

± 

Cb-Co 

^A-Co ; 


(s) 

km/s 

km/s 

km/s 

km/s 

km/s 

km/s ' 

_ 

300.OO 

3.6457 

3.7380 

3.7263 

0.0745 

0.0903 

0.0806 

3 

297.25 

3.7260 

3.7252 

3.7152 

0.0057 

-0.0008 

■0.0108 I 

2 

280.00 

3.6394 

3.6574 

3.6455 

0.0117 

0.0180 

0.0061 

3 

274.73 

3.6350 

3.6367 

3.6244 

0.0057 

0.0017 

-0.0106 

2 

260.00 

3.6213 

3.6013 

3.5872 

0.0613 

-0.0200 

-0.0341 

3 

255.72 

3 5880 

3.5911 

3.5766 

0.0057 

0.0031 

-0,0114 

2 

240.00 

3.5979 

3.5732 

3,5563 

0.0532 

-0.024-7 

-0.0416 

3 

239.30 

3.5690 

3.5724 

3.5556 

0.0058 

0.0034 

-0.0134 

3 

224.90 

3.5680 

3.5692 

3.5503 

0.0064 

0.0012 

-0.0177 

2 

220.00 

3.5611 

3.5712 

3.5516 

0.0352 

0.0101 

-0.0095 

3 

212.12 

3.5790 

3.5746 

3.5538 

0.0071 

-0.0044 

-0.0252 

3 

200.87 

3.5950 

3.5842 

3.5619 

0.0086 

-0.O108 

-0.0331 

2 

200.00 

3.6010 

3.5849 

3.5626 

0.0369 

-0.0161 

-0.0384 

3 

190.34 

3.6150 

3.5957 

3.5723 

0.0U2 

-0.0193 

-0.0427 

2 

190.00 

3.5906 

3.5961 

3.5728 

0.0416 

0.0055 

-0.0178 

3 

180.99 

3.6310 

3,6088 

3,5849 

0.0110 

-0.0222 

-0.0461 

2 

180.00 

3,6365 

3.6102 

3.5863 

0.0232 

-0.0263 

-0.0502 

3 

172.54 

3,6410 

3.6207 

3.5963 

0,0131 

-0.0203 

-0.0447 

2 

170.00 

3.6439 

3.6243 

3.5998 

0.0319 

-0.0196 

-0.0441 

3 

164.83 

3.6520 

3,6316 

3.6067 

0.0168 

-0.0204 

-0,0453 

2 

160.00 

3.6366 

3.6390 

3.6144 

0.0306 

0.0024 

-0.0222 

2 

150.00 

3.6457 

3.6558 

3.6320 

0.0354 

0.0101 

-0.0137 

2 

140.00 

3.6823 

3.6730 

3.6503 

0.0464 

-0.0093 

-0.0320 

2 

130.00 

3.7201 

3.6930 

3.6725 

0.0376 

-0.0271 

-0.0476 

2 

120.00 

3,7329 

3.7150 

3.6976 

0,0303 

-0.0179 

-0.0353 

2 

110.00 

3.7588 

3.7400 

3.7265 

0.0341 

-0.0188 

-0,0323 

1 

110.00 

3 7600 

3.7400 

3.7265 

0.0320 

-0.0200 

-0.0335 

2 

100,00 

3,7105 

3.7694 

3.7603 

0.0470 

0.0589 

0.0498 

1 

100.00 

3.7740 

3.7694 

3.7603 

0.0270 

-0.0046 

-0.0137 

1 

90.00 

3.8120 

3.8040 

3.8000 

0.0220 

-0.0080 

-0.0120 

2 

90.00 

3.7841 

3.8040 

3.8000 

0.0553 

0.0199 

0.0159 

1 

80.00 

3.8580 

3,8433 

3.8442 

0.0230 

-0.0147 

-0.0138 

2 

80.00 

3.7848 

3.8433 

3.8442 

0.0763 

0.0585 

0.0594 

1 

70.00 

3,9000 

3.8837 

3.8896 

0.0190 

-0.0163 

-0.0104 

2 

70.00 

3.8975 

3.8837 

3.8896 

0.0781 

-0.0138 

-0.0079 J 

1 

60.00 

3.9460 

3.9270 

3.9373 

0.0190 

-0.0190 

-0.0087' 

2 

60.00 

3.9860 

3.9270 

3.9373 

0.0920 

-0.0590 

-0.0487 

l!, 

50.00 

3.9980 

3.9703 

3.9831 

0.0230 

-0.0277 

-0.0149 

2 

50.00 

3.9256 

3.9703 

3.9831 

0.1134 

0.0447 

0.0575 

1 

45.00 

4.0200 

3.9914 

4.0046 

0.0230 

-0.0286 

-0.0154 

1 

40.00 

4,0390 

4.0112 

4.0238 

0.0150 

-0.0278 

-0.0152 

1 

36.00 

4.0300 

4.0245 

4.0362 

0.0230 

-0.0055 

0.0062 

1 

32.00 

4.0340 

4.0322 

4.0432 

0.0240 

-0.0018 

0.0092 

1 

28.00 

4.0060 

4.0240 

4,0362 

0.0260 

0,0180 

0.0302 

1 

24.00 

3.9340 

3.9761 

3.9959 

O.OoiO 

0.0421 

0.0619 

1 

20.00 

3.8030 

3.8004 

3.8507 

0.0470 

-0.0026 

0.0477 

1 

18,00 

3.7290 

3.7159 

3.7827 

0.0910 

-0.0131 

0.0537 


References: 

N=1 Mitchell and Yu (i960) 

N=2 Mills and Hales (1978b) 

N=3 Dziewonski and Steim (1982) 
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[ TABLE 5 

j LOVE WAVE PHASE VELOCITIES FOR TWO AVERAGE OCEAN MODELS ’ 
Mi T '■ r ■ l r. i ^ ‘ r r ( r.—r i 



1 

Cn 


C'a 



...Q-Co 


(s) 

km/s 

km/s 

km/s 

km/s 

km/s 

km/s 


300.0 

5.236 

5.230 


0.018 

-0.006 

-0.018 

1 

275.0 

5.157 

5.152 


0.016 

-0.005 

-0.018 1 

1 

250.0 

5.074 

5.074 


0.014 

0.000 

-0.013 i 

1 

225.0 

4.993 

4.999 

4.985 

0.016 

0.006 

-0.008 1 

a 


4.929 

4.926 

4,912 

0.017 

-0.003 

-0.017 

1 

175.0 

4.358 

4.656 

4.841 

0.016 

-0.002 

-0.017 

1 

150.0 

4.791 

4.788 

4.774 

0.015 

-0.003 

-0.018 

1 

125.0 

4.753 

4.725 

4.710 

0.021 

-0.028 

-0.043 

2 

102.4 

4.669 

4.671 

4.655 

0.024 

0.002 

-0.014 

2 

93.1 

4.643 

4.650 

4.634 

0.024 

0.007 

-0,009 

2 

81.9 

4.616 

4.625 

4.609 

0.023 

0.009 

-0.007 

2 

70.6 

4.590 

4.602 

4.585 

0.023 

0.012 

-0.005 

2 

60.0 

4.570 

4.581 

4.564 

0.021 

0.011 

-0.0 0^■ 

2 

51.2 

4.547 

4.564 

4.547 

0.022 

0.017 

0.000 

2 

40.0 

4.527 

4.545 

4.528 

0.023 

0,018 

0.001 

2 

35.3 

4.509 

4.538 

4.521 

0.023 

0.029 

0.012 

2 

30.0 

4.492 

4.530 

4.513 

0.024 

0.038 

0.021 


References: 

N=1 Kanamori (1970) 

N=2 Mitchell and Yu (1980) 
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TABLE 6 

LOVE WAVE GROUP VELOCITIES FOR TWO AVERAGE OCEAN MQDET.S 


4.437 

4.432 

4.412 


11 

15 

20 

25 


4,431 

4.437 

4.444 

4,447 

4.450 

4.453 

4.455 


1 0.032 

0.013 

0.038 

0.021 

0.035 

0.012 

0.030 

0.005 

0.031 

-0.019 

0.034 

-0.009 

n 037 

0.012 

0.U33 

0.031 

0.030 

0.031 

0.029 

0.039 

0.035 

0.061 


Reference: N=2 Mitchell and Yu (1900) 









ORIGINAL PACi% 

OF POOR QUAU'lY 


VELOCITIES FOR 


TABLE 7 
0-20 M.Y. OLD 
RAYLEIGH 


OCEANIC REGION 


km/s 

km/s 

5.239 

5.253 

5.033 

5.049 

4.841 

4.861 

4.676 

4.695 

4.537 

4.552 

4.428 

4.441 

4.335 

4.342 

4,254 

4.255 

4.186 

4.185 

4.123 

4.123 

4.078 

4.076 

4.034 


4,031 

4.025 

3.999 

3.999 

3.999 

3.992 

3.970 

3.970 

3.966 

3.958 

3.944 

3.946 

3.927 

3.927 

3.914 


3.914 

3.913 

3.890 

3.902 

3.900 


3.890 

3.890 

3.875 

3.888 

3.863 

3.884 

3.883 

3.879 

3.886 


3.866 

3.876 

3.868 

3.874 

3.875 

3,875 

3.877 

3.876 

3.877 

3.878 

3.882 

3.881 

3.881 

3.884 

3.879 i 

3.888 


Icm/s 


0.012 

0.011 

0.010 

0.010 

0.009 

0.009 

0.008 

0.008 

0.008 

0.007 

0.007 

0.007 

0.024 

0.007 

0.024 

0.007 

0.025 

0.007 

0.007 

0.024 

0.007 

0,025 

0.007 

0.007 

0.025 

0.009 

0.011 

0.013 

0.025 

0.025 

0.026 

0.027 

0.026 

0.025 

0.025 

0.026 


C-C, I N 


km /s 


0.014 2 

0.016 2 
0,020 2 
0.019 2 

0.015 2 

0.013 2 

0,007 2 

0.001 2 
- 0.002 2 
0.000 2 
- 0.002 2 
- 0.002 2 
-0.006 2 
- 0.001 2 
-0.007 2 

0.000 2 
-0.008 
0.002 
- 0.000 
0.012 
- 0.001 
0.012 
0.001 
0,000 
0.013 
0.001 
-0.004 
- 0.010 
0.010 
0.006 
- 0.000 
- 0.001 
0.001 
- 0.001 
0.003 
0.009 


LOVE 


N 


km/s 


3.602 

3.669 

3.683 

3.683 

3.714 

3,767 

3.809 

3.841 

3.862 

3.884 

3.901 

3.899 

3.889 

3.840 

3.766 

3.512 


km/s 


3.645 

3.666 

3,691 

3.719 

3.747 

3.777 

3.805 

3.820 

3.836 

3.852 

3.872 

3.895 

3.918 

3.931 

3.931 

3.930 


I ± 


< km/s 


0.040 

0.035 

0.028 

0.029 

0.023 

0.024 

0.028 

0.028 

0.027 

0.028 

0.030 

0.033 

0,040 

0.062 

0.118 

0,145 




I km/s 


0.043 

-0.003 

0,008 

0.036 

0,033 

0,010 

-0.004 

- 0,021 

-0.026 

-0.032 

-0.029 

-0.005 

0.029 

0.091 

0.165 

0.418 


300.0 

273.0 

248.0 

225.0 

204.0 

186.0 

169.0 

153.0 

139.0 

126.0 

115.0 

104.0 
102.4 

95.0 

93.1 
86.0 
81,9 

78.0 

71.0 
70.6 

65.0 

60.0 

59.0 

53.0 

51.2 

48.0 

44.0 

40.0 

40.0 

35.3 

30.0 

28.0 
26.0 

24.0 

22.0 
20.0 



km/s 


5.156 

5.056 

5.000 

4.920 

4.845 

4.784 

4.705 

4.641 

4.559 


iz 


km/s 


0.043 

0.051 

0.048 

0.041 

0.040 

0.043 

0.049 

0.044 

0.039 


300.0 

275.0 

250.0 

225.0 

200.0 

175.0 

150.0 

125.0 
102.4 


5.149 

5.069 

4.990 

4.912 

4.838 

4.767 

4.698 

4.633 

4.578 


0.031 

0.031 

0.031 

0.031 

0.031 

0.031 

0.031 

0.031 

0.031 




-0.008 

0.019 


4 374 
4.413 
4.331 
4.363 
4.351 
4.290 


4.335 

4.339 

4.345 

4.350 

4.357 

4.364 

4.372 

4.376 

4.380 
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142.9 

123.0 

111.1 
102.4 
100.0 

93.1 

90.9 

81.9 

77.0 

70.6 

86.7 

60.0 

58.8 

51.2 

50.0 

40.0 

40.0 

35.3 

33.3 

30.0 

28.0 
26.0 

25.0 

24.0 

22.0 
20.0 
20.0 


km/s 


4.254 
4.131 
4.089 
4,065 
4.053 
4.027 
4.025 
I 3.999 
3.998 
3.975 

3.970 

3.959 

3.960 

3,950 

3.952 

3.953 

3.959 
0.958 

3.966 

3.966 

3.969 

3.971 
3.979 

3.971 

3.971 
3.965 

3.959 




0.018 

0.014 

0.019 

0.013 

0.019 

0.010 

0.019 

0.009 

0.019 

0.008 

0.019 

0.008 

0.020 

0,008 

0.020 

0.009 

0.021 

0.021 

0.021 

0.011 

0.020 

0,020 

0.021 

0.012 


I 0.004 
0.007 
0.004 
0.005 
-0.007 
0.002 
0.000 
- 0.000 
- 0.002 
0.000 
- 0.002 
-0.003 
-0.009 
-0.004 
- 0.010 
-0.004 
-0.004 
- 0.002 
-0.009 
0.001 
0.002 
0.005 
0.011 


100.0 

94.3 

90.0 

83.3 

80.0 

73.5 

70.0 
64.9 

60.0 

57.5 

50.8 

50.0 

45.0 

44.6 

40.0 

39.4 

36.0 

34.8 

32.0 

30.8 

28.0 
27.1 

24.0 

23.9 

21.5 

20.0 



TABLE 8 

VELOCITIES FOR 20-50 M.Y. OLD OCEANIC REGION 
RAYLEI GH 


km/s 


4.233 

4.145 

4.097 

4.065 

4.057 

4.034 

4.029 

4.004 

3.991 

3.977 

3.970 

3.959 
3.958 

3.950 

3.950 

3.950 

3.950 

3.954 
3.957 
3.962 
3.965 

3.969 

3.971 

3.971 
3.973 

3.970 

3.970 


I 3.539 
3.684 
3.724 
3.693 
3,714 
3.755 
3.758 
3.776 
3.785 
3.808 
3.828 
3,840 
3.875 
3.884 
3.939 
3.935 
3.969 
3.963 
4.001 
3.986 
4.018 
4.001 
4.022 
3.993 

4.004 

3.977 
3.943 
3.943 
3.893 
3.839 


3.733 

3.750 

3.763 

3.786 

3.797 

3.819 

3.832 

3.852 

3.871 

3.882 

3.910 

3.914 

3.936 

3.938 

3.959 

3.962 

3.978 

3.984 

3.995 

3,999 

4.005 

4.005 

3.991 

3.991 

3.953 

3.901 


km/s 


0.062 

0.045 

0.032 

0.028 

0.027 

0.029 

0,022 

0.019 

0.023 

0.016 

0.019 

0.013 

0.019 

0.012 

0.014 

0.023 

0.023 

0.016 

0.002 

0,018 

0,023 

0.019 

0.024 

0.020 

0.026 

0.020 

0.031 

0.021 

0.025 

0.047 


U-U. 


km/s 


112 
003 
016 
036 
019 
005 
005 
OlO 
0.012 
0.011 
0.004 
0.012 
-0.004 
- 0.002 
-0,029 
- 0.021 
-0,033 
-0.025 
-0.042 
-0.024 
-0,040 
-0.018 
-0.027 
0.006 
0.001 
0.028 
0.048 
0.048 
0.060 
0.062 
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TABLE 9 

VELOCITIES FOR 50-100 M.Y. OLD OCEANIC REGION 
RAYLEIGH 



T 

Q _ 

C 

t 

maoBm 

N 

■BH 




waam 

L is 

km/s 

km/s 

km/s 

km/s 

... 

(s) . 

km/s 

wssim 



fi 

102.4 

4.098 

4,097 

0,012 


1 

110.0 

3.796 

3.766 

0.022 

-0.030 

l 

93,1 

4.075 

4.071 

0,013 

-0,004 

1 

100.0 

3.834 

3.600 

0.019 

-0.028 

1 

81,9 

4.051 

4.051 

0.013 


1 

90.0 

3.066 

3.856 

0.016 

-0.010 

1 

70.6 

4.037 

4.033 

0.013 

-0.004 

1 

80.0 

3.931 

3.907 

0.016 

-0.024 

1 

60.0 

4.028 

4.027 

0,013 


1 

70.0 

3.973 

3.958 

0.013 

-0.015 

1 

51,2 

4.032 

4.029 

0.013 

-0.003 

1 

60.0 

4.016 

4.009 

0.014 

-0.007 

l 

40.0 

4.039 

4.039 

0.014 



50,0 

4.061 

4.051 

0.016 

-0.010 

1 

35.3 

4.044 

4.044 

0.014 

0.000 


45.0 

4,071 

4.067 

0,016 

-0.004 

1 

30.0 

4.045 

4.047 

0.014 

0.002 


40.0 

4.077 

4.077 

0.015 

-0,000 

l 

28.0 

4.045 

4.046 

0.014 

0.001 


36.0 

4.042 

4.077 

0.015 

0.035 

1 

26.0 

4.041 

4,043 

0,014 

0.002 


32.0 

4.047 

4.065 

0.016 

0.018 

1 

24.0 

4.036 

4.035 

0.014 

-0.001 


28.0 

4,007 

4.020 

0.017 

0.021 

1 

22.0 

4.022 

4.020 

0.014 

-0.002 

H 

24.0 

3,924 

3.931 

0.020 

0.007 

l 

20.0 

3.988 

3.990 

0.014 

0.002 

H 

20.0 

3.767 

3,639 

0.030 

-0.128 












d 

T 

c. 

C 

dr 

waam 

N 

T 


u 


mam 


L. (s) _ 

km/s 

km/ s 

km/ s 

km/s 


■Idi 


WS59m 

km/s 

mssim 

1 

175.0 

4.900 

4.901 

0.021 

-0.007 

1 

110.0 

4.496 

4.491 

0.028 

-0.005 

1 

150.0 

4.841 

4.836 

0.021 

-0,005 

1 

100.0 

4.484 

4,496 

0.033 

0.012 

1 

125.0 

4.803 

4.775 

0.021 

-0.029 

1 

90.0 

4.519 

4.502 

0.031 

-0.017 

1 

102.4 

4.743 

4.723 

0.021 

-0.020 

1 

00.0 

4.489 

4.507 

0.026 

0.018 

1 

93.1 

4.716 

4.703 

0.021 

-0.013 

1 

70.0 

4.502 

4.514 

0.026 

0.012 

1 

81.9 

4.694 

4.600 

0.020 

-0.014 

1 

60.0 

4.462 

4.520 

0.028 

0.050 

1 

70.6 

4.669 

4.658 

0.020 

-0,011 

1 

50,0 

4.495 

4.526 

0.031 

0.031 

1 

. 60.0 

4.643 

4.630 

0.010 

-0.005 

1 

45.0 

4.51 

4.529 

0,029 

0.019 

1 

51.2 

4.618 

4.623 

0.019 

0.005 

1 

40.0 

4.499 

4.532 

0,026 

0.033 

1 

40.0 

4.597 

4.605 

0.020 

0.008 

n 

36.0 

4.456 

4.535 

0.024 

0.079 

1 

35.3 

4.584 

4.598 

0.020 

0.014 

D 

32.0 

4,414 

4.536 

0.029 

0.122 

l 

30.0 

4.570 

4.591 

0,021 

0.021 

■1 







Reference: N=1 Mitchell and Yu (1980) 

































OP POOR 


TABLE 10 

VELOCITIES FOR >100 M.Y. OLD OCEANIC REGION 
RAYLEIGH 


102.4- 

93.1 
81.9 
70.6 
60.0 

51.2 

40.0 

35.3 

30.0 

28.0 
26.0 

24.0 

22.0 
20.0 


5.2B8 

5.114 

4.926 

4.760 

4.593 

4,454 

4.374 

4.350 

4.260 

4.146 

4.122 

4.097 

4.085 

4.081 

4,076 

4.072 

4.065 

4.052 

4.042 

4.032 

4.017 

3.986 

3.954 


5.316 

5.130 

4.945 

4.768 

4.606 

4.465 

4.385 

4.339 

4.232 

4.157 

4,133 

4.111 

4.091 

4.081 

4.077 

4.073 

4.070 

4.058 

4.050 

4.038 

4.019 

3.989 

3.935 


0.008 

0.008 

0.006 

0.007 

0.007 

0.007 

0.010 

0.015 

0.040 

0.012 

0.013 

0.013 

0.013 

0.013 

0.013 

0.014 

0.014 

0.014 

0.014 

0.014 

0,013 

0.013 

0.013 


0.009 

0.013 

0,011 

0.011 

- 0.012 

■0.028 

0.011 

0.011 

0.014 

0.006 

- 0.000 

0.001 

0.001 

0.005 

0,006 

0.008 

0.006 

0.002 

0.003 

-0.019 


LOVE 


N 


10.0 

100.0 

90.0 

80.0 

70.0 

60.0 

50.0 

45.0 

40.0 

36.0 

32.0 
2R.0 

24.0 

20.0 


3,953 

4,008 

4.056 

4,071 

4,069 

4,049 

4.018 

3.972 

3.901 

3.781 

3.546 


km/s 


3.843 

3.877 

3.912 

3.950 

3.986 
4.019 

4.040 

4.041 
4,035 
4.019 

3.987 
3.922 
3,773 
3.316 


km/s 


0.013 

0.014 

0.016 

0.015 

0,015 

0.015 

0.015 

0,016 

0.018 

0.026 


km/s 


0.011 

0.006 

0.004 

0.003 

0.022 


-0.038 

-0.031 

-0.028 

-0.015 

0.001 

0.015 

0.021 

-0.008 

-0.230 


4,586 

4.560 

4.539 

4.510 


4.512 
4,531 
4.517 

4.478 4.461 0.022 

4.490 4.467 0.023 

4.504 4.473 0.025 

4.477 4.478 0.027 

4.448 4.481 0.024 

4.424 4.484 0.022 

4.441 4.485 0.021 

4.438 4.486 0.025 


Re ferences: 

N=1 Mitchell and Yu (1980) 

N=2 Souriau and Souriau (1983) 
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TABLE 11 

VELOCITIES FOR 0*5 M.Y. OLD 
RAYLEIGH 


OCEANIC REGION 


4.338 

4.131 

4.072 

4.012 

3.967 

3.925 

3.005 

3.043 

3.022 

3.003 

3.794 

3.784 

3.780 

3.771 


4.190 

4.1ls6 

4.056 

4.009 

3.974 

3.921 

3.085 

3.859 

3.833 

3.807 

3.794 

3.786 

3.789 

3,792 


0.007 

0.041 

0.028 

0.025 

0.022 

0.017 

0.014 

0.012 

0.012 

0.013 

0.015 

0.016 

0.021 

0.027 


C- 


km/s I 


0.140 

-0.015 

-0.017 

-0.003 

0.007 

-0.004 

0.000 

0.016 

0.011 

0.004 

0.000 

0,002 

0.001 

0.021 


fs) 

km/s 

km/s 

km/s 

■ km/s 

137.0 

3.476 ! 

3.605 , 

0.095 

i 0.129 

122.0 

3.584 ' 

3,615 

0.069 

1 0.031 

107.5 

3,643 

3.624 

0.041 

-0.019 

94.3 

3,613 

3.632 

0.041 

0.019 

83.3 

3.625 

3,642 

0.028 

0.017 

73.5 

3.643 

3.651 

0.023 

0.008 

64.9 

3.650 

3,662 

0.019 

0.012 

57.5 

3.668 

3.672 

0.018 

0.004 

50.8 

3.694 

3.683 

0.020 

-0.011 

44.6 

3.729 

3.695 

0.023 

-0.034 

39.4 

3.753 

3.700 

0.026 

-0.045 

34.8 

3.756 

3,724 

0,027 

-0.032 

30.0 

3.762 

3.742 

0.028 

-0.020 

27.1 

3.769 

3.762 

0.028 

-0.007 

23.9 

3.774 

3.782 

0.031 

0.008 

21.5 

3,829 

3.797 

0.040 

-0.032 

10.7 

3.762 

3,005 

0.047 

0.043 

16.4 

3.828 

3,013 

0,077 

-0.015 


Reference: Forsyth (1975a) 
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TABLE 12 

VELOCITIES FOR 5-10 M.Y, OLD OCEANIC REGION 
RAYLEIGH 


N 

T 



i: 

c-c„ 


^"■"T — 


'"'X'' ' 


U_~^LV: 


fs) 

km/s 

km/s 

km/s 

km/s 


(s) 

km/s 

km/s 

km/s 

km/s 1 


142.8 

4.174 


0.100 

0.026 

1 

137.0 

3.493 


0.132 

0.107 ’ 

1 

125.0 

4.044 

4.117 

0.062 

0.073 

1 

122.0 

3.568 

3.615 

0.108 

0.047 i 

1 

111.1 

3.999 

4.056 

0.043 

0.057 

1 

107.5 

3.753 

3.629 


-0,124 ; 

1 

100.0 

3.961 

4.011 

0.041 

0.050 

1 

94.3 

3.672 

3.646 

0.075 

-0.026 i 

1 

90.9 

3.975 

3.976 

0.038 


1 

83.3 

3.648 

3,666 

0.047 

0,020 

1 

77.0 

3.930 

3.927 

0.028 


1 

73.5 

3.636 

3.685 


0,049 

1 

66.7 

3.890 

3.896 

0.025 


1 

64.9 

3.690 

3.707 


0.017 

i 

58,8 

3.862 

3.875 

0.021 

0.013 

1 

57.5 

3.723 

3.730 


0.007 

1 

50.0 

3.846 

3.856 

0.021 

0.010 

1 

50.8 

3.793 

3.754 


-0.039 1 

1 

40.0 

3.848 

3.844 

0.022 


1 

44.6 

3.810 

3.782 


-0.027 

1 

33.3 

3.852 

3.845 

0.025 


1 

39.4 

3.827 

3.814 


-0.013 

1 

25.0 

3.886 

3,863 

0.031 

-0.023 

1 

34.8 


3.849 

0.049 

-0.036 1 

1 

20.0 

3.903 

3.887 

0.036 

-0.017 

1 

30.8 

3.992 

3.886 

0.051 

-0.106 

1 

16.7 

3.888 

3.898 

0.043 

0.010 

1 

27.1 

3,925 

3.925 

0,052 

0.000 1 







I 

23.9 

3.875 

3.967 

0.054 

0.086 







1 

21.5 

3.752 

3.983 


0.231 







1 

18.7 

3.694 

3.993 

0,074 

0.299 







1 

16.4 

3.830 

4.005 


0.175 1 


Reference: Forsyth (1975a) 
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1 

! Uppermantle velocities 

TABLE 15 
For each Mltchel 

! thicknesses 

* 0-80 

1 80-50 

^ A 

1 water 

! 3.45 

1 4.67 

5.05 

sediment 

0.02 


0.18 

crust 1 

1.51 

' 1.58 

1.60 

crust 8 

4.64 

5.15 

5.19 

lid 1 

8.00 

8.00 

6.00 

lid 8 

! 14.00 

24.00 

! 89.00 


WaSM 

■BBEa 



km/s 

km/s 

km/s 

water 

1.52 

1.52 


sediment 

1.65 

1.65 

1.00 

LID 1 

8.08 

8.19 

4.40 

crust 1 
0-80, 80-50 

5.21 

5.21 

3.03 

B 

5.15 

5.15 

3.00 

A 

4,94 

4.94 

2.88 

50-100 

5.07 

5.07 

2,96 

>100 

5.01 

5.01 

2.93 

crust 8 
0-80, 80-50, B 

6.80 

6.80 

3.90 

A 

6.53 

6.53 

3.74 

50-100 

6.70 

6.70 

3.84 

>100 

6.83 

6.63 

3.80 

LID 8 
0-20 

8.81 

6.21 

4.80 

80-50 

8.42 

8.42 

4.72 

A, B 

8.41 

8.41 

4.71 

50-100 

8.39 

8,48 

4.70 

>100 

8.27 

8.31 

4.63 

LVZ top 
0-80 

7.67 

7.90 

4.20 

80-50 

7.77 

7.88 

4.28 

B 

7.48 

8.00 

4.36 

A 

7.92 

8.02 

4.36 

50-100 

8.04 

8.15 

4.43 

>100 

8.12 

8.12 

4.48 

LVZ bottom 
0-80 

7.57 

7.77 

4,31 

20-50 

7.59 

7.77 

4.32 

B 

7.53 

8.05 

4.43 

A 

7.52 

7.87 

4.28 

50-100 

7,46 

8.03 

4.25 

>100 

7.66 

8.03 

4.36 

880 

0-20, 20-50 

8.47 

8.47 

4.80 

A. B, 50-100. >100 

8.56 

8.56 

4,64 

400 

0-20, 20-50 

8.82 

8.82 

4.72 

A. B. 50-100, >100 

8.91 

8.91 

4.77 
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TABLE 18 


Comparison of two average o cean models 
"1 Hi ^ <150 I >150 


! SW 

S 

SW 1 s 

SW 

S 

Model 1 
ua 
Ur 

11.55 

26,92 

11.40 

18.51 

8.238 

33,98 

13.90 

22.43 

2.449 

35.51 

13.51 

18.51 

2.630 

35.37 

16.87 

22.43 

12.44 

17.64 

3.63 

8.762 

27.86 

3.86 

Model 2 
Or 
% 
ol 

Ur. 

10,70 

21.58 

11.95 

_13.49_ 

7.257 

26.73 

14.99 

22.42 

6,340 

26.99 

15.71 

18.49 

6.327 

28.94 

17.51 

22.42 

10.47 

15.18 

3.681 

7,449 

24,65 

4.010 

Totals 
Model 1 
Model 2 

24.15 

18.79 

24.72 

21.29 

28.10 

25.37 

20.50 

24.54 

14.75 

10.99 

19.01 

16.55 



FIGURE CAPTIONS 


Figure. 1: 


Figure. 2: 


Figure. 3: 


Figure. 4: 


Figure. 5: 


Phase velocity dispersion curves for the average Earth, Lower curve 
is for Rayleigh waves, Upper curve is for Love waves. Data points are 
indicated by +. Where no error bars are visible standard deviations 
are Less than or equal to the width of the horizontf=»l bar of the +. 
Solid lines give the best fitting average Earth model. 

Group velocity dispersion curves for the average Earth. Symbols as in 
figure 1 

Phase and group velocity dispersion curves for the average ocean. 
Lower curve in each plot is for Rayleigh waves, upper curve Is for 
Love waves. Where no error bars are visible the standard deviations 
are about half the width of the point for phase velocity, and the order 
of the line thickness for group velocity. The solid lines give the 
dispersion calculated for average ocean model A. 

Same as figure 3 but the solid lines give the dispersion calculated for 
average ocean model B. 

Phase and group velocity dispersion curves for Mitchell and Yu’s 0-20 
M.Y. old oceanic province are given in the leftmost two plots. 
Corresponding plots for their 20-50 M.Y. old oceanic province are 
given in the rightmost two plots. Deta..s of are the same as for figure 
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Figure.6; Phase and group velocity dispersion for Mitchell and Yu's 50-100 M./ 
old oceanic province are given in the leftmost two plots. 
Corresponding plots for their >100 M.Y. old oceanic province are 
given in the rightmost two plots. Details are the same as for figure 3. 

Figure.7: Phase and group velocity dispersion for Forsyth's 0-5 M.Y. old oceanic 

province are given in the rightmost pair of plots. Corresponding plots 
for 5-10 M.Y, old oceanic province are given in the leftmost pair of 
plots. Details are the ssme as for figure 3. 

Figure. 8: Phase and group velocity dispersion for Forsyth's 10-20 M.Y. old 

oceanic province are given in the rightmost pair of plots. 
Corresponding plots for his 0-10 M.Y. old oceanic age province are 
given in the leftmost pair of plots. Detaiiu of each plot are the same 
as for figure 3. 

Figure. 9: Changes in Rayleigh wave phase velocity dispersion curves caused by 

changes in each of the parameters varied in the study. Except where 
otherwise noted changes are in the LVZ. From right to left the 
parameters changed are the following: TOP ROW: SH velocity 
gradient(SKG), mean SH velocity (SHM), ETA gradient (ETAG), mean 
ETA (ETAM), lid thickness (THL). MIDDLE ROW: SH and PH velocity in 
the lid with SH/PH constant (LH), mean PH velocity (PHM), PV 
velocity gradient (PVG), PH velocity gradient (PHG), SV velocity 
gradient (SVG). BOTTOM ROW: mean PV velocity (PVM), crustal 
velocity (VCR), crustal thickness (THCR), SV and PV velocity in the lid 
with SV/PV constant (LV), mean SV velocity(SVM). Horizontal scale is 
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frotn Os Lo 300s 

Figure, 10: Changes to Rayleigh wave group velocity dispersion curves due to 
changes in each of the parameters varied in this study, The 

parameters are displayed in the same order as figure 9, 

Figure. 11; Changes in Love wave phase velocity dispersion curves due to 

changes in each of the parameters varied in this study. Except where 
otherwise indicated changes are in the LVZ. 

From right to left the parameters changes are the following; TOP 
R0W:PVM, PVG. ETAM, ETAG, PHG. MIDDLE ROW; LV, SVM, SHG, SVG, 
THL. BOTTOM ROW: SHM, VCR, THCR, LH, PHM. 

Codes used are those in the caption of figure 9, Horizontal scale is 
Os - 300s. 

Figure. 12: Changes in Love wave group velocity dispersion curves due to 

changes in each of the parameters varied in this study. The 

parameters are displayed in the same order as figure 11. 

Figure, 13: Velocity depth profiles for two average ocean models that fit the data 
equally well. Dispersion curves for model A (solid lines) are shown in 
figure 3. Dispersion curves for model B (dashed lines) are shown in 
figure 4. From left to right the top row gives SV, PV, and ETA and the 
bottom row gives SH, and PH. 
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Figure, 14: Rayleigh wave phase and group velocity residuals for the two average 
ocean models. Upper right plot shows model B group velocity 
residuals, lower right shows model B phase velocity residuals. Upper 
left plot shows model A group velocity residuals and lower left shows 
model A phase velocity residuals. A positive residual indicates that 
the model is too fast. 


Figure. 15: Love wave phase and group velocity residuals for the two average 
ocean models. Layout is the same as figure 14. 

Figure. 16: Velocity depth profiles for the average Earth model. From left to 
right solid lines show PV, SV, and ETA, and dotted lines show PH, and 
SH. Dispersion curves for this model are shown in figures 1 and 2. 

Figure, 17: Velocity depth profiles for the 0-20 M.Y, (upper set) and the 20-50 
U,Y. (lower set) old age provinces of Mitchell and Yu. Details of each 
set are the same as figure 16. Corresponding dispersion curves are 
given in figure 5. 

Figure. 18: Velocity depth profiles for the 50-100 M.Y, (upper set) and the >100 
M.Y. (lower set) old oceanic regions of Mitchell and Yu. Details of 
each set are the same as figure 16. Corresponding dispersion curves 
are shown in figure 6, 
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Figure. 19: Velocity depth profiles for the 0-5 M.Y, (upper set) and the 5-10 M.Y. 

(lower set) old oceanic regions of Forsyth, Details of each set are the 
same as in figure 16, Corresponding dispersion curves are shown in 
figure 7. 

Figure. 20: Velocity depth profiles for the 10-20 M.Y. (upper set) and the 0-10 
M.Y. (lower set) old oceanic regions of Forsyth. Details of each set 
are the same as in figure 16. Corresponding dispersion curves are 
shown in figure 8. 

Figure, 21: The thickness of the seismic lithosphere as determined in this study. 

The upper edge of the open boxes gives the thickness of the LID only. 
The lower edge gives the thickness of the LID plus the crust. The 
elastic thicknesses and the isotherms are from Watts et. al. [1900]. 
The triangle is a refraction measurement of lithosphere thickness 
from Nagumo et. al. [1981]. 

Figure, 22: Schematic representation of seismic velocities due to temperature, 
pressure, and crystal orientation assuming a flow aligned olivine 
model. The upper left diagram shows a convection cell with arrows 
indicating flow direction. The trench Is indicated by T, the ridge by 
R, and the midpoint by M. The lower left diagram shows temperature 
depth profiles for the trench, ridge, and midpoint. The upper and 
lower right diagrams show the nature of the velocity depth structure 
of VSH and VSV respectively due to pressure, temperature, and 
crystal orientation. 
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